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Abstract 

We study nonperiodic tilings of the line obtained by a projection method with an inter- 
val projection structure. We obtain a geometric characterisation of all interval projection 
tilings that admit substitution rules and describe the set of substitution rules for each 
such a tiling. We show that each substitution tiling admits a countably infinite number of 
nonequivalcnt substitution rules. We also provide a complete description of all tilings of 
the line and half line with an interval projection structure that are fixed by a substitution 
rule. Finally, we discuss how our results relate to renormalization properties of interval 
exchange transformations (with two or three intervals). 
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1 Introduction 



The study of substitution rules for symbolic sequences and one-dimensional tilings has a long 
and rich history. In particular, Sturmian sequences and their substitution rules have been 
intensively studied in the context of computer science and dynamical systems [Pyt02, Lot02]. 

Hedlund and Morse [MH40] showed that Sturmian sequences can be given an elegant geo- 
metric construction (see also [Ser86, Ser85]): every Sturmian sequence can be obtained from 
the intersections of a straight line with a square grid by assigning the symbol to each intersec- 
tion of the line with a horizontal grid line and a 1 to each intersection with a vertical grid line. 
The binary cutting sequence thus obtained is a Sturmian sequence. Sturmian substitution rules 
are those morphisms of binary strings that map Sturmian sequences to Sturmian sequences. 
These substitution rules have an elegant geometric realization: a unimodular matrix with inte- 
ger coefficients inducing a linear transformation of the line associated with the initial Sturmian 
sequence to the line associated with its image under the substitution rule. Interestingly, the 
Sturmian substitution rules are generated by only three morphisms. Incidentally these three 
morphisms also generate the automorphism group of the free group on two symbols F2, see eg 
[PWWOO, Lam98, WW94]. 

In the physics literature, in the area of quasicrystals [Sen95, AG95], Sturmian sequences 
arise as symbolic representations of canonical projection tilings. Their construction, which we 
describe in Section 2, is closely related to the one by Hedlund and Morse mentioned above. 

In this paper we consider one-dimensional interval projection tilings. This class of tilings 
appears to have been first considered in [MPPOOb], and contains the above mentioned set of 
canonical projection tilings. 

Our results unify and extend various results in the literature. In particular, we answer the 
following central questions: 

• Which one-dimensional tilings with an interval projection structure are substitution 
tilings? 

• What substitution rules do such tilings admit? 

We note that the first problem was previously only resolved for one-dimensional canonical 
projection tilings. The second question had not been addressed even in the case of particular 
well-studied examples such as the Fibonacci tiling. 

Whereas many studies of nonperiodic tilings have been mainly combinatorial or algebraic, 
in this paper we take a predominantly geometric point of view. 

Importantly, it turns out that the results of this paper provide some key insights enabling 
the understanding of the existence of substition rules for higher dimensional (canonical) pro- 
jection tilings, including the well-known Penrose and Ammann tilings [Har03, HL04] which 
serve as a prototypical examples of quasicrystalline tilings. 

2 Main results 

We consider tilings of the line by intervals obtained by the following geometrical construction. 
Let V and W be two transversally intersecting lines in M^. We consider the intersection of the 
lattice C with a strip F -|- C M^, where O is a half open interval of the vertical axis, 
referred to as the window. For later convenience, we adopt the convention that the window 
lies in the subspace Y parallel to the direction the vertical lattice generator. Subsequently, we 
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project the lattice points inside Z"^ n {V + to V, choosing the projection Hy parallel to W 
(such that Uy^ {v) is a Hne parallel to W for all v eV). The resulting discrete point set 

Uv{z'^n{v + n))cv 

is considered as the set of vertices for a tiling of the line V, where tiles are intervals. The 
construction is sketched in Figure 1. We refer to the resulting tilings as one- dimensional 
tilings with an interval projection structure. In this construction we always assume that V and 
W are not parallel to one of the coordinate axes (thus avoiding some degenerate constructions 
yielding periodic tilings). 




Figure 1: Illustration of the geometric structure of the Fibonacci tiling, which is the prototypical 

example of a onc-dimcnsional tiling with an interval projection structure. Consider a strip in 
R^, indicated in grey and defined by the translation of a unit square in the direction of a line 
V with slope t~^, where r = (1 + -\/5)/2 denotes the golden ratio. The intersection of with 
this strip is projected down the space ly to y to yield the set of vertices of the Fibonacci tiling 
on V. The window is the interval of the space Y that generates the strip when translated along 
V. 

In the special case where the length of the window fl is precisely equal to the length 
of the projection Ily to Y (parallel to F) of a 1 x 1 unit square in R^, the points inside 

the strip can be connected by horizontal and vertical line pieces (steps) to form a monotone 
staircase. This staircase can be viewed as a discrete approximation of a line parallel to V. 
Such staircases can be represented by sequences of two symbols, each symbol representing a 
type of step (horizontal or vertical). Such binary sequences are the Sturmian sequences and 
have been extensively studied in the literature [Lot02], also in relation to substitution rules 
[dB81, Sec98, IY90. Par97, Ber96, CMPS93]. Sturmian sequences arise in number theory as 
the sequence of differences between consecutive terms in a Beatty sequences (the numbers 
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{[{nX\))^^i where A is irrational). For discussions of substitution rules in this context see 
[Sto76, FMT78, dB89, KvdP96, Kom99]. 

When W is chosen such that the projections parallel to W onto V of consecutive steps of 
the staircase intersect only at their endpoints, the tiling of V can be represented by a binary 
sequence; each symbol representing projection of either a horizontal or a vertical step. 

In general, we consider tilings of the line by intervals (of finite length), where the intervals 
are taken from a finite protoset of tiles (sometimes called prototiles). Every tile in the protoset 
has a geometric shape (an interval with a given length) and possibly a label. A tiling consists 
of an infinite collection of translated copies of tiles from the protoset that cover the line in 
such a way that translated prototiles only overlap on their boundaries. Labels provide a way 
to distinguish between geometrically identical tiles in a protoset. 

In this paper we provide a characterisation of all one-dimensional tilings with an inter- 
val projection structure that admit a substitution rule. We also characterise the set of all 
substitution rules that such a tiling admits. 

The (geometric) substitution rules that we consider consist of two parts. First, we perform 
a uniform inflation of each (labeled) tile by a constant factor, say A > 1. Thereafter, we replace 
each inflated tile by a set of uninflated tiles in such a way that each inflated tile with the same 
label is replaced in exactly the same way. Here, we always require that the substitution is 
vertex hierarchic, implying in the present context that each inflated tile is exactly covered by 
the set of tiles that replace it. We say that a tiling is a substitution tiling if the image of 
a tiling under a substitution rule is locally isomorphic to the original tiling. Two tilings are 
locally isomorphic if they are on any local scale indistinguishable in the sense that any finite 
patch of the one also appears in the other. A set of locally isomorphic tilings is called a local 
isomorphism class. In the case of one-dimensional tilings with an interval projection structure, 
local isomorphism classes are characterized by the spaces, V and W and the window length 

(see Lemma 3.6). 

A well known example of a one-dimensional substitution tiling with an interval projection 
structure is a Fibonacci tiling. Its construction is as sketched in Figure 1 with V having its 
slope equal to the inverse of the golden ratio (r~^), W = V-^, and the window CI having length 
1 + T~^. The corresponding tiling can be represented by a binary Sturmian sequence and 
admits the well known Fibonacci substitution rule. A geometric derivation of the Fibonacci 
substitution rule is presented in Figure 2. 

In this paper we establish how the geometric approach to a substitution rule for the Fi- 
bonacci tiling, as illustrated in Figure 2, generalizes to provide a full characterization of all 
substitution rules for one-dimensional substitution tilings with an interval projection structure. 

It should be noted that we view the tilings in the first instance as geometric objects (a 
partition of the line by intervals) without considering labels. Such a geometric tiling admits 
a substitution rule if there exists a labelling of the tiles (respecting the tile-shapes, ie two 
intervals of different length cannot be assigned the same label), such that the labelled tiling 
admits a substitution rule^ . 

^In particular this rules out certain Sturmian sequences that are substitutive with respect to symbolic sub- 
stitution rules (ie they can be generated by projection on the symbols). This is because the projection breaks 
the geometric requirement that tiles of different lengths cannot have the same label. For example consider the 
substitution rule a — » ac, 6 — » acacb, c — » acadb. When this substitution rule is brought into the geometric 
setting the tiles come in two lengths, the a tiles and the h and c tiles. As the substitution rule does not dif- 
ferentiate between h and c, we can remove the labels to obtain a Sturmian sequence that can be constructed 
using a geometric substitution rule. However if we bring together the a's and the 6's we get a different Sturmian 
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Figure 2: Geometric derivation of a substitution rule for the Fibonacci tiling. 

(a) A patch of the Fibonacci tihng and the corresponding staircase. The tihng space V and 
W = (along which we project to are the expanding and contracting eigenspaces of the 



matrix M = 




(b) The vertices of the staircase that have a vertical step to their right are identified. They 
correspond to lattice points lying in the dark grey strip. 

(c) After removal of these points, we obtain a staircase with horizontal and diagonal steps. We 
may interpret this removal as the substitution a6 — * a, a ^ 6, where a represents a long and h 
represents a short tile. We note that in the orginal tiling horizontal staircase steps project to 
long tiles and vertical steps to short tiles, whereas in the new staircase horizontal steps project 
to short tiles and diagonal steps to long tiles. 

(d) By application of the lattice automorphism M~^, the new staircase is mapped to a staircase 
for a Fibonacci tiling (as characterized by the slope of V and the same window length), as it 
commutes with the projections along V and along W . It follows that the tiling for the new 
staircase is a Fibonacci tiling that is uniform^ inflated by a factor equal to the expanding 
eigenvalue of M. The substitution rule a ^ a6, 6 ^ a is known as the Fibonacci substitution 



The first result of this paper is summarized in the following theorem: 

Theorem 2.1. Let V and W be the tiling and window spaces for a nonperiodic one- dimensional 
tiling with an interval projection structure. This projection tiling admits a substitution rule if 

and only if V and W are the expanding and contracting eigenspaces of a primitive matrix in 
G/(2,Z), and the window length is inside Q[A] where X is the expanding eigenvalue of M ? 

It is important to note that the above theorem deals with nonperiodic tilings. Periodic 
tilings are special, as these may admit substitution rules that bare no relationship with the 
geometric setting of Theorem 2.1. For instance, if a one-dimensional tiling with an interval 
projection structure is periodic with unit cell U, and W is chosen such that the ratios of the 
lengths of the tiles are rational, then there exists an appropriate inflation (multiplication by 
some integer) such that all inflated tiles can be replaced by a number of concatenated unit 
cells U. 

In Theorem 2.1 no reference is made to the protoset. In fact, it turns out that each 
one-dimensional substitution tiling with an interval projection structure admits a countable 
infinity of substitution rules, using many different labelings of the tiles (and thus different 
protosets). The following theorem provides a geometrical characterization of all substitution 
rules of one-dimensional substitution tilings with an interval projection structure. 

Theorem 2.2. Consider a nonperiodic one- dimensional substitution tiling with interval pro- 
jection structure and window fi, as in Theorem 2.1. Then for each choice of the (half open) 
interval window ^' (ZW such that 

(i) = 

(ii) M{Q + V) C + V), 

(Hi) the length of each connected component of Q,' \ U.yM{Q,) is in Q[X], 

there exists a substitution rule that maps the projection tiling with window to the locally 
isomorphic projection tiling with window il' . 

Moreover, for each substitution rule one can identify a window Cl' with the above mentioned 
properties. 

Among the substitutions mentioned in Theorem 2.2, we may distinguish between local and 
nonlocal substitutions. Local substitutions are those where the shape of a finite neighbourhood 
of a tile (a patch consisting of a finite number of adjacent tiles) determines its label. In 
many studies of substitution rules, such as [PleOO, BSJ91], the locality of substitution rules 
is assumed. However, the following result shows that although one-dimensional substitution 
tilings with an interval projection structure always admit infinitely many nonlocal substitution 
rules, only some admit local substitutions. 

sequence. This symbolic projection, however, gives the same label to tiles of different lengths, thus we do not 
consider this to be a geometric substitution rule. For a general discussion of substitutive Sturmians from a 
symbolic point of view see [AR91]. 

^Recall that Q[A] = Q + AQ is the smallest algebraic field containing A, Gl{2, Z) is the group of 2 x 2 matrices 
with integer coefficients whose determinant is equal to ±1. A matrix M € Gl{2,Z) is primitive if their exists 
n e N such that all entries of M" are non-zero. 
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Theorem 2.3. Consider the characterization of substitution rules in Theorem 2.1 and The- 
orem 2.2. Let the window pair and Q' be associated to a substitution rule for a tiling with 
interval projection structure. Let M G Gl{2,Z) be the primitive matrix associated to this sub- 
stitution rule. Moreover, let d+ and d- denote the distances between the tops and bottoms, 
respectively, of the windows nyM(r2) and Q! . 

Then the substitution rule is local if and only if at least one of the following conditions is 
satisfied: 

• {det{M)-X)\n\,d+,d-eZ[X], 

• (det(M)-A)|Q|,Ad+-det(M)|0|,Ad_-det(M)|fi| eZ[A]. 

Corollary 2.4. A substitution tiling with interval projection structure whose window length 
satisfies 

' ' ^ detM + A detM- A^ 
only admits nonlocal substitution rules. 

To the best of our knowledge, the existence of nonlocal substitution rules was not known 
before, not even for the widely studied Fibonacci tilings, for which we provide an example. 

Example 2.5 (Nonlocal substitution rule for Fibonacci tilings). Let the tiling space 
V and projection direction W be the expanding and contracting eigenspaces of the matrix 

M-{\ ,1, 

The Fibonacci tiling is obtained by the projection method, with a window that is obtained as 
the projection of a unit square along V to Y. 

We consider the substitution rule where the expansion predecessor is associated with the 
subwindow M{fl) that lies precisely at the centre of the initial window Q' . It is readily verified 
that the distance between the endpoints of the subwindow M{Q) and the endpoints of W is 
T~^/2 G Q[A] \ (Z[A]). Hence, by Theorem 2.3 the corresponding substitution rule is nonlocal. 
It is defined on a protoset consisting of four tiles: two long intervals (ai and 02) and two short 
intervals (61 and 62)- It terms of these, the substitution rule takes the form 




The derivation of this substitution rule is presented in Example 7.2. 

Example 2.6 (Substitution tiling with only nonlocal substitution rules). Consider a 
projection tiling with V and W the expanding and contracting eigenspaces of the matrix M 
of (1), as in the case of the Fibonacci tiling. We consider tilings with window length 1/n with 
n > 1 integer. We have A = -r and detM = -1 and 1/n Z[r]/(1 + r) U Z[r]/(1 - r), so 
that by Theorem 2.1 and Corollary 2.4 the tiling admits substitution rules but none of them 
is a local substitution rule. 
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Historically, there has been a considerable interest in tilings that are fixed by a substitution 
in the sense that the tiling, after the application of a substitution rule, is an exact (translated) 
copy of the original tiling.^ 

The following theorem summarizes our results about the existence of one-dimensional sub- 
stitution tilings with an interval projection structure that are fixed by a substitution rule. 

Theorem 2.7. Consider the set of substitution rules S fixing the local isomorphism class 

of a one- dimensional substitution tiling with an interval projection structure associated to a 
window length \Q\ and matrix M € G/(2,Z). Then, a substitution rule in S fixes det(M — /), 
det(M — /) — 1, or dct(M — /) + 1 tilings, where I denotes the identity matrix. 

The conditions on a substitution rule that determine which one of the three options given 
in the above theorem arises, can be expressed in terms of properties of the windows fl, CI' and 
the matrix M. We refer the reader to Theorem 6.1. In Section 6 we also describe the set of 
one-dimensional tilings with an interval projection structure that tile the half-line M"*" and are 
left invariant by a substitution rule. 

Our geometric methods are inspired by and adapted from the algebraic methods used by 
Masakova et al. in [MPPOOb]. Masakova et al. present a large class of quasicrystal tilings of 
the line. Their tilings have an interval projection structure where V and W are the expanding 
and contracting eigenspaces of a primitive matrix M with integer coefficients. Masakova et al. 
use a different, but equivalent, description of these tilings in terms of algebraic number theory. 
Under the assumption that M is unimodular they show that such a tiling admits a labelling and 
a substitution rule on the labelled protoset that fixes the tiling, if and only if both endpoints 
of the window lie in Q[A] in Y. This can be viewed as a corollary of the results in this paper. 

Our results can also be viewed in the context of renormalization properties of two and 
three-interval exchange transformations. For a discussion, see Section 7. 

As the proofs of our results are all constructive, they can be viewed as the basis of an 
algorithm for constructing substitution rules. In Table 1 we have listed some examples of 
substitution rules found using the methods described in this paper, encoded as a Mathematica 
program [HL]. 

It turns out that the substitution rules in Table 1 (and all others we obtained) are invert- 
iblc, in the sense that it has a formal inverse so that it corresponds to an automorphism of a 
free group. It is well known that binary Sturmian substitution rules correspond to automor- 
phisms of the free group F2 [WW94, PWWOO, Lam98, EI98]. Our experiments thus suggest 
that all substitution rules for one-dimensional tilings with an interval projection structure are 
invertible. 

3 Preliminaries 

In this section we recall some definitions and properties of substitution rules and interval 
projection tilings. 

3.1 Substitution rules 

Consider a set of prototiles r = {tj}, and a real number A > 1. Let At = {Arj} denote the 
set of tiles obtained from r after inflating each prototile by a uniform factor A. A (vertex 

^Wc adopt the convention that two tilings are equivalent if they can be mapped onto each other by a 
translation. 
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M 



S2 



Substitution 



Inverse 



1 1 

2 1 



a biab2 
bi bia 
b2 ab2 



bi 
b2 



62 a" 
ab. 



2 

-1, 



1 1 

2 1 



ai — 6i6iai 

02 ^ 6103^2 

03 ^ &1 0-3^1 

62 —>■ 61 ai 



02 ^ 

03 ^ 62ar^"3^20.r^ 



6201 ^bi 



61 ^ Olbn 



0162 ^^3 ^02 



a ^ca 



1 1 

2 1 



a ^ be 
b ^ ba 
c — > 6cc 



-1, 



b ^ ac "a 
c 



-1, 



1 1 

2 1 



a 



c 



-^cb 

cbab 



be ^ba2 ^cb ^ 
bc-^b 



2 3 
1 1 



ai 60102 
02 —>■ ca2 
b — ^ caia2 
c — caicaia2 



Oi 



02 



b aib-^a2b-^dr^ 
c — 026^^06"^ 



Table 1: Some examples of substitution rules for one-dimensional tilings with an interval pro- 
jection structure, where V and W are the expanding and contracting eigenspaces of a matrix 
M G Gl{2, Z), and the width of the window O is represented by I € Q^: |0| = IIy{1). The sub- 
window associated with the substitution rule is represented by the vectors si,S2 € Q^: IIy{si) 
and IIy{s2) are the relative positions of the bottom and top of the subwindow. Note that 
Ml = Sl — 32- The labels for the tiles are chosen such that those labelled corresponds to the 
projection of a horizontal staircase step, those labelled bi correspond to the the projection of a 
vertical staircase step and those labelled q correspond to the projection of a diagonal staircase 
step (the combination of a horizontal and vertical step). In the last column we provide the 
formal inverse of the substitution rule. 



hierarchic) replacement rule ur associates a patch of tiles to each expanded prototile in Ar, 
so that aii{XTi) covers the inflated prototile Atj exactly. Consequently, the replacemt rule is 
vertex hierarchic: the vertices satisfy the inclusion vert{Ti) C verta^^Xri). 

A substitution rule a consists of a pair (X,aR), indicating the subsequent application of 
an uniform inflation by the factor A and a replacement rule aR. The action of a substitution 
rule (7 on a protoset r induces, a substitution rule on any tiling (or patch) T admitted by this 
protoset. A patch or tiling is substituted by applying the uniform inflation, giving AT, and 
replacing all translated copies of the inflated prototiles according to the replacement rule. 

A tiling T' is called a predecessor (tiling) for T if a (7"') = T. We refer to the inflated tiling 
XT' as the expansion predecessor (tiling). 

Definition 3.1 (Substitution Tiling). A tiling T, with protoset r, is a substitution tiling 
with substitution rule a if it has inflnite predecessors under a, ie if for all n G N there exists a 
locally isomorphic tiling T' such that cr^T' = T. 

It follows that the set of substitution tilings for a primitive^ substitution rule forms a set 
''a substitution rule is called primitive if after a finite number of subsequent applications each tile (in the 
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of local isomorphism classes. 

The geometric definition of substitution tiling given here is motivated by the study of 
geometric substitution rules in higher dimensions, see [Sen95, GS98]. In one dimension these 
are closely related to symbolic substitution rules and morphisms as studied in computer science 
and dynamical systems, see [AS03, Pyt02]. 

The process of finding a substitution rule for a tiling T is a two stage process: identifying 
candidates for the expansion predecessor and then examining which of these can be related 
to T by a replacement rule. An expansion predecessor tiling T must be locally isomorphic to 
AT, where A is the expansion factor of the substitution rule. Moreover, as we consider vertex 
hierarchic substitution rules, any candidate for expansion predecessor must also satisfy the 
inclusion vert{T) C vert{T). 

Having found a tiling T that satisfies these two conditions, the question is whether there 
exists a replacement rule relating T to T, or equivalently if each tile in T is covered by the same 
set of tiles in T. The latter step possibly requires assigining labels to distinguish between tiles 
with the same geometric shape (interval of a given length) . Finally, it must be verified that the 
substitution rule on the labelled protoset admits a substitution tiling (infinite predecessors can 
be found). It is important to note that this last fact does not directly follow from the existence 
of a substitution rule relating T to T: it needs to be shown that T has a predecessor with 
respect to the same substitution rule, and more generally the existence of an infinite sequence 
of predecessors. 

It may well happen that one identifies a substitution formulated on too large a protoset, 
due to unnecessarily many labels for certain geometric tile shapes. This observation leads us 
to discuss a natural notion of equivalence between substitution rules. We propose here the 
notion of patch equivalence, under which two substitution rules are taken to be equivalent if 
for any number of iterations of the substitution rule on the tiles from the protoset, the sets of 
obtained patches (while ignoring the labels) are the same. In order to formalize this notion, let 
vr denote the function that extracts the geometric shape of a labelled prototile (interval with 
given length) or labelled patch (sequence of adjacent intervals). 

Definition 3.2 (Patch equivalent substitution rules). Let r and r' be two labelled 
protosets and a and a' be substitution rules on respectively r and r'. Then a and a' are 
patch equivalent if for every prototile T e t and all n G N there exists a T' G r' such that 
7r(cj"^T) = 7r(((j')"^T'), and similarly if for every prototile T' G r' and all n G N there exists a 
T G r such that 7r(a"r) = 7r((a')T'). 

We call a substitution rule patch minimal if there exists no patch equivalent substitution 
rule with fewer prototiles. It is readily verified that within an entire set of patch equivalent 
substitution rules for a substitution tiling, there exists a unique patch minimal one. 

protoset) is replaced by a set of tiles including copies of all tiles in the protoset. As a consequence, under the 
subsequent application of a primitive substitution rule no patches grow up that contain only a subset of tiles 
from the protoset. 

®Note that in order to render this local isomorphism class unique, other (less intuitive) definitions of substi- 
tution tilings arc used in the literature, sec for example [GS98]. It turns out that the projection structure used 
in this paper forces the substitution rules to be associated to a unique local isomorphism class. 
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3.2 One-dimensional substitution tilings with an interval projection struc- 
ture 

Consider a one-dimensional tiling with an interval projection structure. The latter means that 
that the tiled line can be embedded in M? in such a way that the vertices of the tiling are 
projections from points of the lattice I? that lie within a bounded distance of the embedded 
line. These lattice points may be linked together by line segments in to form a staircase, 
as illustrated in Figure 3. 




Figure 3: Illustration of a tiling with interval projection structure and its associated interval 
projection staircase. The tiling is constructed by the projection to the subspace V of a. slice 
of the lattice Z^. Here V here is the line with gradient j-^^ (the inverse of the golden ratio) 
and the slice is defined by the intersection of the lattice with a strip parallel to V . The lattice 
points within this strip are projected orthogonal, along the subspace W = . By connecting 
lattice points within the strip whose projections yield neighbouring vertices in the projection 
tiling by line segments, we obtain the interval projection staircase. The window is the subset 
of the space Y that gives the slice when translated along V . 

The vertices of this staircase can be projected to the window space Y that is transversal 
to V . The neighbourhood of a vertex within the staircase (and the neighbourhood of the 
corresponding vertex in the tiling) can be deduced entirely from the position of its projection 
along V to Y . For instance, one can partition the window in regions indicating which staircase 
step lies to the right of a vertex whose projection to Y lies in any of these regions. In Figure 4, 
one finds such a window partition for the three possible right steps of the staircase in Figure 3. 

The example in Figure 3 displays three different staircase steps, projecting to three different 
interval tiles in V . It turns out that this number of staircase steps (and tile shapes) is in fact 
typical for one-dimensional tilings with an interval projection structure. This is the content 
of the next Lemma, which also features in Masakova et al. [MPPOOa, MPPOOb] and is related 
to the three distance theorems in number theory. For a historical account of the latter, see 
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Strip 



Strip for a 



Strip 




Figure 4: The window partition (and associated partition of the strip) for the three possible 
staircase steps (to the right) of the interval projection staircase in Figure 3. These three types 
of steps give rise to three geometric tile types in the projection tiling. Note that the boundaries 
of the strip partition are obtained by translating the edges of the strip by two of the three 
staircase steps. The corrsponding window partition is obtained by projection to Y. 

[AS03, pp53-54]. We include a short and self-contained proof for completeness. 

Lemma 3.3. A one-dimensional projection tiling with an interval projection structure has tiles 
of at most three lengths. In the case of three lengths, one is equal to the sum of the two others. 

Proof. Consider the points Hyit) and Ily (t + I) at the ends of the interval window fli^t C W 
with length ITy (/), where I G and base f G M^, and the translated lattice + 1 (which has 
a point in common with the base of ^i^t)- Now consider the set of points {^i^t + ^) n {I? + t), 
and its projection (along W) to V. 

Since Cli^t is bounded, this projection yields a discrete point set in V. We now consider the 
two points in I? + t that project to either side of ny(t) in V. We denote these points t + a 
and t — b, where a,b correspond to steps of the staircase. 

For simplicity, we continue the proof assuming that the window is the half-open interval 
^l,t = \t^t + I) d W. The argument for the remaining case in which ^ = {t,t + l\ d W \s 
analogous. 

The subwindows corresponding to these two steps are the intervals fi"^ C Q,i^t and ^\^C 0,i^t 
given by 

Qft = {peni,t\'nY{a)+penu} = [TlY{t),UY{t + l-a)), 

nl = {peni^t\^Yib)+peni^t} = [^Yit-b),UY{t + i)). 

If ny(t + I — a) = Hyit — b) then all steps in the staircase are accounted for, and there are 
only two types of steps in the staircase. 

If ny(i + I — a) ^ ny(i — 6), the situation is as shown in Figure 5. Apart from the 
subwindows for a and b discussed above, there is a remaining third subwindow in between. 
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Y 




ny(t + o 



ny(a) ^ 

Figure 5: The partition of the window for a staircase with three steps. The arrows indicate 
how nearest neighbour staircase steps project to the window. The central subwindow is the 
subwindow for staircase step a + b. 

Consider a point g G whose projection to Y lies in this central subwindow, i.e. ny(q) G 
[Ily (t + I — a), Ily (t — b)). Then the projection to Y of the point q + a + b lies in the interval 
[Ily (i + / + 6), Ily + a)) which is a subinterval of fli^f We finally claim that if g G projects 
to the central subwindow, the next (right) step in the staircase is a + 6. Namely, if there are 
points whose projection to V lies between Ily{q) and IIv{q + a + b), then by periodicity of the 
lattice t + a and t — b would not project to the closest points to t, contradicting the definition 
of a and b. 

Hence, there are three types of staircase steps: a, b and a + b. The corresponding sub- 
windows are nf^ = [ny(t),ny(t + I - a)), ^1+^ = [Uyit + I - a),ny(t - b)) and Of^ = 

[UY{t-b),UY{i + l)). 

In the case of the other choice of the half-open window interval Q^t = (ny(t),ny(i + I)], 
one finds nf^ = (ny(t),ny(t + I - a)], 9."+^ = (ny(t + Z - a),ny(t - 6)] and = (ny(t - 

6),ny(t + z)]. ' ' □ 

The staircase can be seen as the orbit of a function ^i t^ which maps a vertex of the staircase 
to its right nearest neighbour. 

Definition 3.4 (The staircase function The function ^i^t '■ {V + ^l,t) {V + Q^) 

is defined as 

' z + a if Uriz) G ^If^, 



z + a + b if ny(z) G 
z + b ifny(z)G^^'' 



where Q^j, f^"^'', and are the subwindows of O;^^ for the steps a, a + b and b, as introduced 
in the proof of Lemma 3.3 above. 

As the function ^i^t is invertible, by iterating one can also generate successions of 
neighbours to the left, ^i^t is defined for every point r inside the strip V + 0;^^, so that the set 
{$"j(r) I n G Z} contains the vertices of an interval projection staircase in r -|- Z^. 

This is just a translation of the interval projection staircase in Z^ with window ^i^t-Uyir)- 

The subwindows for a, b and a + b can be further subdivided, providing a partition relating 
projections of vertices of the staircase to patches (collections of adjacent tiles) instead of single 
tiles (or staircase steps). 

Let us illustrate this in more detail. 
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Consider the interval [Ily (t), Ily (t+Z— a)) in Y. The vertices that project to this subwindow 
are followed by the step a. Call the corresponding set of vertices of the staircase Sa- Because 
of the definition of Sa we have ^i^tiSa) = Sa + a, representing the set of vertices at the end of 
an a step. The projection Uyi^i^tiSa)) densely fills the window interval [Ily (t + a), Ily (t + 1)), 
the subwindow for a translated by ny(a). 

For each of the vertices in Sa + a, the next right step (after a) depends on the subwindow 
to which it projects. If one or both of the points ny(i + I — a) and ny(t — b), lie in [ny(i + 
a), Ily (t + /)) then the latter interval can be partitioned into several subintervals, on the basis 
of their projection to one of the three subwindows for a, b and a + b. The image of these 
subintervals under yields a partition of the subwindow for a, signalling the next two steps 
that follow a vertex (the first one being a). 

This process can be continued by induction to construct window partitions, corresponding 
to sequences of steps of the staircase (and corresponding tiles in the projection tiling) after a 
given vertex. The points in the window that act as boundary points for these partition intervals 
are the images of the two boundary points under iterates of the map ^^l'- Ily ($j^"(ny (t+/ — a)) 
and Ily ($[""(ny (i— 6)). We may therefore associate tiling (and staircase) patches with window 
intervals. 

Definition 3.5 (Window Interval for the Patch P). The window interval Ip for the patch 
P of a onc-dimcnsional tiling with interval projection structure is the unique interval of the 
window such that IIy{z) £ Ip for some vertex z of the projection staircase, if and only if the 
vertex 11^(2;) of the tiling is followed (to the right) by the patch P. 

We now use the window intervals to show that the set of one-dimensional tilings with an 
interval projection structure with fixed V, W and I forms a local isomorphism class. 

Lemma 3.6. The set of one- dimensional interval projection tilings characterised by the linear 
subspaces V and W ofM? satisfying = VF, and fixed I G M^, forms a local isomorphism 
class of one- dimensional tilings. 

Proof. Let T^t denote the interval projection tiling with window vector I and base t. We first 
show that the tilings Ti t smd Ti t' ^ire locally isomorphic for any t,t' G M^. Consider a patch 
P of Ti^ti with window interval Ip C O^^j. The equivalent window interval Ily (/p — t + t') can 
be found in fi; ^'. Any vertex of the staircase for TJ ^/ which projects to this interval will be 
followed by the patch P, and as ny(Z^) is dense, such vertices exist. Thus every patch of Ti t 
occurs in Ti^t' , by an analogous argument any patch of 7/ occurs in Ti^t so the two tilings are 
locally isomorphic. 

We now show that any tiling of the local isomorphism class has an interval projection 
structTire. Consider a tiling Tthat is locally isomorphic to Ti t- The tiling must have the same 
tile shapes so we can find a staircase S, with the same steps as the staircase for T^t, such that 
11^(5) = T. Now consider the closure of the projection of S to Y: Ily (5) C W. If the extreme 
points of this set were more than ny(/) apart, then there would be two vertices of S laying 
more than ny(Z) apart in the projection to Y. This is impossible, as these points would lie 
at the ends of some finite patch, which would appear in Ti^f We need only show, therefore, 
that the set Ily (5) is an interval. To do this consider the projection of longer and longer finite 
lengths of the staircase to Y. These have smaller and smaller gaps between the points as each 
finite length of staircase is also part of the staircase for Ti t, whose window is an interval. We 
may therefore get arbitrarily close to any point in the interval between the extreme points of 
Ily (5) so that the closure yields the interval. 
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This concludes the proof of the fact that all tilings in {Ti^t\t G M^} are locally isomorphic 
and any tiling locally isomorphic to Ti^t is contained in the set {^,t|i € M?}. 

□ 

Note that if ny(t) = Hyit' + z) for some z G then Ti t and 7/,^/ differ only by a 
translation. Thus, with respect to the natural equivalence relation (see footnote 3) each value 
of t G y/ny(Z^) is associated with a unique tiling. 

4 The partition algorithm 

In this section we discuss substitution rules in the context of the geometry of the interval 
projection structure. Section 4.1 provides a discussion of the setting, after which in Section 4.2, 
we describe the existence of substitution rules, in this context. In Section 5 we examine 
precisely when substitutions rules exist, and provide a chactcrisation of all substitution rules 
for any given tiling of the line with an interval projection structure. 

The proofs of our results are all constructive, and in Section 4.3 we illustrate the construc- 
tion of a substitution rule in the case of a specific tiling. The algorithm has been implemented 
in Mathematica and is available from the authors [HL]. 

4.1 Window partitions and replacement rules 

In this section we show that every one dimensional tiling with a interval projection structure 
and substitution rule can be associated to a unimodular 2x2 primitive integer matrix. The 
action of this matrix on the projection structure is key to providing a geometric interpretation 
of the substitution rule. 

We first associate a matrix to any substitution tiling with an interval projection structure. 
This procedure is folklore, see for instance [Ken90]. 

Lemma 4.1. The spaces V and W of a nonperiodic^ one- dimensional projection tiling with 
interval projection structure with a substitution rule are the eigenspaces of a non-singular 
primitive, 2x2 integer matrix. 

Proof. Consider a tiling T with interval projection structure from the lattice C and substitution 
rule a. Let S be the staircase for T. Now consider the expansion predecesssor T for T under 
a. Wc may lift T to a staircase S whose vertices form a subset of the vertices of S. This 
staircase lies on a lattice JC which is a sublattice of C The steps in this lattice corresponding 
to the tiles a and 6 in 7" generate this lattice, just as the equivalent steps in S generate C We 
may therefore define the linear map taking £ to £ as the map taking the steps corresponding 
to tiles a and 6 in 5 to the equivalent steps in S. 

The action of this matrix must commute with the projection, thus the spaces V and W are 
the eigenspaces of the matrix. As the tiling is non-periodic, the space V must be irrational, 
and the matrix must therefore be primitive and non-singular. □ 

A key result of this paper concerns the fact that the matrix mentioned in the previous 
lemma must be unimodular, and thus represents an automorphism of the lattice in the embed- 
ding space. In previous studies [PleOO, MPPOOb, LGJJ93] unimodularity of the matrix was a 
hypothesis. 

^Recall that nonperiodic refers to the geometry of the tihng before adding any labelhng. 
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A substitution rule consists of two components: inflation and replacement. First the tiling 
is homogeneously inflated by a constant factor. In the embedding space we may represent this 
inflation as multiplication by a matrix M, as identified in Lemma 4.1. This maps the staircase 
S for the tiling T to MS which serves as a staircase for the inflated tiling. 

In the embedding space the replacement rule corresponds to replacing every staircase step 
in MS by a section of staircase steps with vertices on the lattice C. 

We now discuss how the substitution rule acts on the window. 

The existence of a substitution rule for an interval projection tiling implies the existence of 
a patch minimal substitution rule for the same tiling. So we assume without loss of generality 
that (J is patch minimal: for some iteration a"^ of a the prototiles can be distinguished by the 
patches they substitute to. Each of these patches is associated to a particular subwindow of 
CI (where Q is the window for T). 

Let T be the expansion predecessor of T with respect to substitution rule cr". Every 
instance of a prototilc T in T is replaced by the same patch P of tiles in T. However P can 
also arrise at other places in the tiling where it is not covered by T. For example P could 
occur as a subset of a patch replacing a different prototile. This occurs in the Fibonacci tiling 
with the substitution rule (a ab,b ^ a): the tile b in the expansion predecessor is replaced 
by the tile a, but a is also contained in the patch ah that replaces a. Thus the subwindow flp 
for P within Q is not necessarily the subwindow $7^ within $7 associated to the prototile T, 
but it must contain it. 

By Lemma 3.3 the window Q may be divided up into two or three subwindows corresponding 
to the lengths of the tiles. Let ^\t\ be the subwindow for all tiles of the same length as T. The 
intersection Q^^l H $lp corresponds to every tile of the same length as T in T that is covered 
by P. Consequently, we have ^It = ^\t\ H Op. 

The replacement rule takes the instances of T and replaces each of them with the same 
patch of tiles. In the staircase the corresponding steps in S are each replaced by a collection of 
steps {si, . . . , Sn} in S. The vertices of the Sj's project to Y (after closure) as the subintervals 
nr(si) + of the window O. Furthermore, as a'^ is a substitution rule, each ny(si) + Oy is 
contained within the interval associated to a unique prototile of T. 

Theorem 4.2. The spaces V and W of a nonperiodic one- dimensional projection tiling with 
interval projection structure with a substitution rule are the eigenspaces of a primitive matrix 
in Gl{2,Z). 

Proof. From Lemma 4.1 we have a matrix M G gl{2,'L), it remains to be shown that M is 
unimodular. 

The vertices of the staircase for the tiling T are the points in which project to the 
window O in y. The vertices of the staircase for the expansion predecessor T are the vertices 
of Ml? which project to the subwindow O. The lattice MI? is a subgroup of 1? and both 
are Abelian groups. We may therefore consider 1? as a finite union of distinct translations of 

Ml?, 

1? = MI? U {gi + MI?) U . . . U {gN-i + MZ^). 

We denote the corresponding set of translations as G = {0, 51, . . . , ^jv-i}. Note that N = 
|det(M)|. 

The replacement rule acting on T yields T. By the discussion above we may view the 

replacement procedure in the projection to the window VI <Z Y . The vertices of a labelled 
predecessor staircase step project to the points IIy{MI?) n [a,b] C O. The replacement rule 
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induces translations of this interval by elements of ny(Z^). Consider one such translation, by 
IIy{z), for some z in Z^. We may associate gi £ G to z where z G MJ? + Qi. This translation 
adds to the vertices of the staircase the points in Ml? + Qi that project to [a, 6] + PiY{z) 
in Y. In order to get every point in Ylyi^'^) that projects to $7, therefore, the action of the 
replacement rule on the subwindow must cover the window at least N times. Furthermore 
the set of intervals covering any region must include an interval translated by an element of 
JIy{MI? + gi) for each gi G G. 

Let V be the set of intervals induced by the action of replacement rule on the labelling 
intervals of the predecessor. For g e G, let the set of intervals in V translated hy z e 
I{y{AH? + g) be Vg. The labelling partition of the window is defined by a finite number 
of points {vi, . . . , Vp}. The subwindow Q, is partitioned by the points {tJi, . . . ,Vp\ in the same 
manner. Each translate of a labelling interval in the subwindow partition induced by the 
replacement rule corresponds to a unique labelled tile and is thus the subset of an interval of 
the partition on Jl. Consider vi. As the end point of an interval of the partition on $7, it must 
be the end point of any interval in V that contains it. Consider the intervals for which v\ is 
a left end point. These must contain an element of Vg for each g E G. The left end points of 
these intervals are the translations of points in IIy{{vi, . . . ,Vp}). Thus for each g E G there 
exists Vg £ {vi, . . . ,Vp} and Zg £ Ml? such that vi = Vg + ny(zg + g). Furthermore, Vg is 
equal to if and only if g = h, as elements of G are not in MZ^. 

The subwindow associated to the q-th iterate of the substitution rule is partitioned in the 
same manner as O by p points. The absolute value of the determinant of the matrix for this 
substitution rule is N*^. Thus the lattice is partitioned into N'^ translates of M'^l'^. Thus if 

> 1, there exists q such that A^'^ > p. This is leads to a contradiction, as Vg is equal to 
Vh if and only \i g = h, so every element of g must have a distinct Vg. Hence N = 1 and 
consequently M is unimodular. □ 

The following lemma describes the set of potential expansion predecessors: tilings, whose 
vertices are a subset of the original tiling, and locally isomorphic up to scaling. 

Lemma 4.3. Let M be a primitive matrix in Gl{2, Z), with eigenspaces V and W and expand- 
ing eigenvalue A. Let il/ j be the window of the interval projection tiling Ti t with window ami 
tiling spaces V and W and t G . Let Qi^t,s = ilj^o + fly i'^ + t), be the window of the interval 
projection tiling Ti^t^g, where s G such that C O/^j. The interval projection tilings XTi i 
and Ti^t,s (^fe locally isomorphic for all t G M^. Furthermore vert{'T i^t,s) C vert{Ti^t)- 

Proof. By Lemma 3.6 the tilings {Ti^t\t £ M^} are locally isomorphic. 

We show that every tiling in {'Ti^t,s\'t £'^,s £ M^} is contained in {\Ti,t\t € Let Si^t 
be the staircase associated to Ti t. As the action of M commutes with the projection, we have 
Jlv{MSi,t) = XTi^t and IW{MS^) = iA-^i^,,*. Hence 

U^,s = HyMO,,o + Hy(s) + Hy(t) = HyMO^^n^ (.)+t 

Thus T^,s = yTi,s+t, so every tiling in {Ti^t,s\t£ M^, s £ M^} is a tiling in {XT^\t G K^}- 

Finally, if ^i^t,s C 0;,^ then every vertex of ^ ^ is a vertex of Ti f. □ 

4.2 The algorithm 

Lemma 4.3 describes the set of potential expansion predecessor tilings. It remains to show for 
which of these there exists a replacement rule. Consider the window Q with subwindow Q, for 
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a potential expansion predecessor. Both fl and i7 are partitioned into two or three intervals 
associated to the two or three tile types given by Lemma 3.3. The tiles of T are covered by 
patches of T. These are finite patches of T so they have interval windows. The subwindow 
n is therefore partitioned into intervals each of which is associated to a specific tile type in 
the expansion predecessor with a specific patch of tiles that cover it. We use this partition 
to label the tiles in T. The identification of an appropriate rescaling of 0, with induces a 
corresponding labelling partition of O. In turn this provides a labeling of the tiles in T. 

In the procedure above we initially identify tile types only by length. We now repeat the 
procedure, initially identifying tile types by length and label. This produces a labelling of the 
tiles, that is either the same as the previous labelling or a refinement of it. In the former case 
we have identified a labelling with respect to which the tiling admits a substitution rule. In 
the latter case we may repeat the procedure. The tiling admits a substitution rule for this 
potential expansion predecessor if and only if after a finite number of iterations no further 
refinement of labelling occurs. 

To describe this procedure in detail, we reintroduce the points I, t and s that describe the 
position of the window fli^t and the subwindow 0,i^t,s- We define a function gi^t,s '■ V + ^l,t 
V + nit by: 

where 

n = min{n G Z+jHy (^^-"(z)) G n^,s}- 

and is the staircase function defined in Definition 3.4. 

To illustrate the action of this function, consider x & 7? CiV + Qi^f Using the staircase 
function we follow points on the staircase until we reach a point p that projects to the 
subwindow ^i,t,s- We apply an affine transformation, based on M~^, that takes V + ^i^t,s 
to V + This affine transformation relates points in Jl; j + V and ^i^t,s + whose Ily- 
projections have the same relative position in their respective windows. The action of the 
function gi^t^g is illustrated in Figure 6. 

Consider the points t + 1 — a and t — b, whose Ily-projections give the boundary points of 
the partition of the window fi; t associated to tile type, as in Lemma 3.3. Let Gi^t,s denote the 
subset of Cli^t to which all gi^t,s iterates project, 

Gi,t,s = ^Y{{gr,t,sit + I - a)\n G Z+}) U UYi{gr,t,sit - b)\n G N}). 

As, Gi^t,s = Gifi^s + ny(t), the relative position of Gi^t,s as a subset of 0;^^ is independant 
of t. " ' ' 

Finiteness of Gi^t,s is a necessary and sufficient condition for the existence of a substitution 
rule. 

Lemma 4.4. The set Gi^t,s is finite if and only if there exists a substitution rule for the interval 
projection tiling Ti^t, with expansion predecessor T; ^ ^ associated to the subwindow C 0,^- 

Proof. We first prove that there is a substitution rule if the set Gi^t,s is finite. 

If Gi^t,s is finite, it provides the boundary points of a partition of fij^t into a finite number 
of intervals. Consider such an interval (51, 52)- In the subwindow ^i^t,s the equivalent interval 

is {gT,g^) = M{{gi,g2)~t) + UY{s + t). 

We show that (51 , ^) n Gi,t,s is empty. Note that 

9l,t,s{i9i,92)) = (91,52) 
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(a) (b) 



Figure 6: Illustration of the action of the function gi^t,s to a point x G V + Cli^f (a) We consider 
a point p in the V + Qi.i and then apply to walk the point back until the image lies in 

^i,t,s + V. Here and do not lie fli,t,s + V but does, (b) Subsequent 

application of the affine transformation. 



as (51,52) C ni^t,s- Therefore 

Since {{gl^g2) n Gi^t,s) is empty by construction, it follows that (51,52) H Gi^t.s is empty. 

As ny(t + Z - a) and nY(t - 6) are in Gi^t^s-, and (5r,52) n Gi^t,s is empty, (5r,52) must lie 
entirely in the subwindow associated to one of the tile types defined in Lemma 3.3. Denote 
the associated staircase step c ^ 1? . Consider the translation of (51,52) by ny(c) (ny(5i + 
c),ny(52 + c)). We show that (Ily (51+ c), Ily (52 + c)) n G;^^^^ is empty. Namely, suppose that 
this set contains a point p. Under gi^t,s^ P is taken to p — c by which is in the subwindow 

^l^t,s-i and then mapped to M~^{p — c—t) + s+t in $7; ^. This yields a point in {{gl, g2)r]Gi^t,s)j 
which is empty by construction. 

By induction of the above argument, there exists an such that for all positive n < N, 
UYO<^f^{{gl,g^))nGi^t,s is empty, and IlYO<^^^{{gl,g^)) C ^^,3- This implies that the partition 
of Cli^t with boundary points in Gi^t,s is the set of labelling intervals for a substitution rule. 

Wc now prove the converse: if Gi^t,s is infinite there is no substitution rule. 

Suppose there is a substitution rule a, with an associated labelling partition of the window 
Ul^f Let the vertices of this partition be the set of points P. The set P must contain the points 
Ily + 1 — a) and Ily — b). The set Gi^t,s is the union of the gi^t,s orbits of Ily + 1 — a) and 
Ily (t — b), one of which must be infinite. As P is finite there must exist some point p £ P such 
that gi,t,s{p) ^ P- Consider the labelling interval that contains gi,t,s{p)- The corresponding 
interval I in the subwindow ^m,s is translated to a series of intervals of the window $7; i, by the 
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replacement rule for a (as discussed in more detail above). One of these translated intervals, 
contains p in its interior. As p € P, it is on the boundary between two labelling intervals. 
Thus we can find two tiles whose label is associated to / that are covered by two different 
labelled patches after applying the replacement rule. This contradicts the assumption that a 
is a substitution rule. □ 



4.3 Illustration 

In this section we illustrate the partition algorithm for finding a labelling partition for a 
substitution tiling and then contruct the corresponding substitution rule. 

We consider the tiling with interval projection structure introduced in Figures 3, 4 and 6 
above. In this example: 




with V and W the expanding and contracting eigenspaces of M. We have, ny((l,0)"^) = r/r/, 
ny((0, 1)^) = 1/r], ny((l,0)^) = -1 and ny((0, 1)^) = r. Recall that r is the golden ratio 
and define 77 = y/2 + r. The window has length (2 + 2r) in Y and the staircase has steps 
a = (1, —1)"^, b = (—1, 2)'^ and a + 5 = (0, l)"^. From Figure 4 one observes that the window is 
divided into three intervals with boundary points Uyit + l — b) = ny(t + (—1, 0)^) = ny(t) + 1 
and livit - a) = Hy (i + (-1, 1)^) = Hy (t) + (1 - r). 

Our aim is to describe the construction of the substitution rule associated to the subwindow 
'^i,t,s = [ny(t) + l,ny(t) + 3) = liviMVlifi + t + s) with s = (-1,0)'^ of the window Vli^t = 
[t, t+{2 + 2r)), with I = (-2, 2f . See also Figure 6. As the substitution rule is independent 
of t, we set t = (0, 0) without loss of generality. 



4.3.1 The partition algorithm 

The algorithm is illustrated in Figure 7. 

We start with the points —a and I — b, which project to W as the boundary points of the 
partition for tile type, —a = (—1, 1)'^ and I — b = (—1,0)"^. We calculate gi^t,s{—a)- the point 
ny(— a) is in Qi^t,s so a) = M^^{—a — s) = (—1,2)-^. Next, we calculate gf^^{—a). We 

have ny((-l,2)^) i Il^,s- Moreover <^-^\i-l,2f) = (0,0)^ and ny((0,0)^) ^ n^,s, but 
2)-^) = (—1, 1)-^ = —a projects to ^i^t,s- Hence —a = gf^g{—a) and the gi^^t^s-orbit of 
—a has period two. 

We repeat the procedure for / — b. We have ny(/ — 6) G 0,i^t,s- Thus gffg{l — b) = 
M~^{1 — b — s) = (0,0)^. Subsequently gf^gil — 1,2)^ = gi,t,s{—0'): considered above. 
Thus the tail of the (7i,t,s-orbit oil — h coincides with the gi^t^s-oihit of —a. 

The union of the ^^^t^g-orbits discussed above contains four points whose Hy-projection 
yields Gi^t,s- In Figure 7 these points are labelled "Window points". By applying M and then 
translating by s, the set of window points are mapped to the set of points labelled "Subwindow 
points" in the figure. The relative positions of the Hy-projection of the subwindow points in 
^l^t,s is in correspondence with the relative positions of the points Gi^t,s in ^i,t- 

The points Gi^t.s induces a partition of $1/ j into four labelling intervals. One interval is 
associated to the tiles of length Ilv{b), one to the tiles of length Ilv{a + b), and two to the 
tiles of length Ilv{a). We refer to the corresponding labelled prototiles as ai, 02, a + 6 and b. 
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M-^-a-s) =5o,s(-a) 
M-i($-2(M-H-a - .)) - .) = gU-a) 



-a, $-2(M-i(-a-.)) 




s, I — b 

(o,of,$,;/(M-H-a-.)) 





Subwindow points 


« 


Window points 




Window /Subwindow 




Other noints 





Figure 7: Illustration of the construction discussed in Section 4.3.1. The arrows indicate the 
action of 

4.3.2 The replacement rule 

By following the translations of the labelling intervals of fli^t,s induced by we obtain the 
replacement rule for the tiling and subwindow considered above. To illustrate this we consider 
the replacement rule for the tile a + 6 in Figure 8, showing that a + 6 is replaced by (0+6)01602- 
For details we refer to the figure. 

In a similar fashion one can obtain the replacement rules for the other tiles, yielding the 
substitution rule: (oi — (o + 6)01,02 — 01602, 6 — ^ (o + 6), o + 6 — ^ (o + 6)01602). 



5 Substitution rules 

In this section we prove our main results about the characterisation of interval projection 
substitution tilings and their substitution rules. 

The following lemma describes when the set Gi^t,s (that has been defined in Section 4.2) 
is finite. Recall from Lemma 4.4 that this condition is central to the characterisation of 
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ITy of this hne 
is the window _ 
for a + 6 steps in 
the subwindow 




Figure 8: Finding the replacement rule for a + b. The interpretation of strips and points is as in 
Figure 7. The figures illustrate the start of the $;^t-orbit of the subwindow of fli^t,s associated to 
the tile a + b. Wc note that every iterate of this interval lies entirely within the strip assocatcd 
to one of the labelled tiles (ai, a2, a + b or b), until the third iterate, in (d), returns to ^i,t,s- 
The $;,t-orbit visits the strips corresponding to the tiles (a + b)aiba2 in order. It follows that 
a + 6 is replaced by (a + b)aiba2. 

substitution rules. 

Lemma 5.1. Gm,s is finite if and only Ily (Z), Ily (s) G Q(A). 

Proof. We first show that Gi^t,s is finite if ny(/),ny(s) G ny(Q2). Note that ny(Q2) = Q[A]. 
Since only the projections to y of Z and s are meaningful parameters for the tiling, we may 
assume without loss of generality that l,s £ Q^. We show that the set {(7q ^(Ily — a))|n G N} 
is finite. By a similar argument it can be shown that the set {g^ g{IlY{t — 6))|n G N} is finite 
as well. Recall that the set Gi^t,s = Gi^^g + t, so from here we only consider Gi,o,s, without loss 
of generality. 

As each staircase step has positive length and a + 6 is the largest possible step, we have 



Now consider the value of n in the definition of g^t,s- As the steps in the expansion tiling have 
a maximum length of ny(M(a + 6)) and the steps in the original tiling have a minimum length 
of min{ny (a), Ily (6)}, n has an upper bound 



Uv{z -a + b)< ny($-/(z)) < Uv{z). 



It thus follows that 



Uv{z-r{a + b)) < ny($-; 



{z))<Uv{z) (reN). 




Now let 7j = glQg{IlY{t + 1 — a)). Prom the observations above, and the definition of $;,o 
and gifi^s we obtain 

Uv{ji-i - N{a + b)) < nv{M-fi + s)< Uv{-fi-i), 
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since Hy (M7i + s) = Ilv{Mgifi^s{n) + s) = Ily ^(7i_i)) and n< N. Consequently, as M 
commutes with the projection Ily, we have 

ny(M-i(7,_i - N{a + 6) - s)) < Uv{^i) < Uv{M-\^i_i - s)). 

By induction on the inequalities above we obtain the following inequality 

i 

ny(7i) > Uv{M-\jo) - {J2 M-')(iV(a + 6) + s)). 

k=0 

The matrix acts on V as multiplication by < < 1, so that for all i 

ny(7o - + + ^)) ^ ^vili) < ny(7o). 

Now recall that s G Q^, thus s e Z^/rUs for some mg G N. We observe that for any n G N 

gifiAi^V^sn) n (y + o,,o)) c ((zVm.n) n (y + o,,o))- 

As 7o = G we thus have 

ji G ZVm, n {no + [Uvijo - + + s)),Uv{jo)]) 

for all i Thus 7^ is contained within a finite subset of the lattice 1? /rng- As I G Q^, we have 
/ G 1? /mi for some m; G N. 

Using an analogous argument one can show that the points generated by applying gi^^s to 
I are also finite. The set G;^o,s is therefore the projection of the union of two finite sets and is 
therefore finite. 

It remains to be shown that whenever Gi^o,s is finite we must have ny(Z),ny(s) G Q[A]. 
We have gifi^l) G M~^{1 — s) + 1?. As W is an eigenspace of M, we have 

ny(5lo,.(0)e vny(z-s) + ny(z2). 

If the set G;,o,s is finite there must exist j G N with i ^ j, such that 

- X^)UYil - s) G ny(Z2). 

However, from the properties of M we have A G ny(Z^) and it follows that if G;^o,s is finite we 
must have ny(^ — s) G Q[A]. The same argument can be repeated, starting with the action of 
5;^o,s on (the other extremum of the window), yielding ny(s) G Q[A]. In turn it follows that 
also ny(Z) G Q[A] □ 

Our main result is the following Theorem, from which Theorem 2.1 and Theorem 2.2 are 
obtained as immediate corollaries. 

Theorem 5.2. Let V and W be the tiling and window spaces for the set of interval projection 

tilings {Ti^t\t £ I^^}; with fixed window vector I (and window length IIy{1))- These interval 
projection tilings are interval substitution tilings if and only if V and W are the eigenspaces 
of a primitive matrix M G Gl{2,Z), and IIy{1) G Q[A]. 

Each substitution rule is associated to a subwindow J^;,t,s = IlyMJl^^o + ny(s + t) with 
s G such that ni^t,s C Substitution rules associated to different choices of s are not 

patch equivalent. 
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Proof. By Lemma 4.4 Gi^t,s is finite if and only if there exists a substitution rule, and by 
Lemma 5.1 the set Gi^t,s is finite if and only if the conditions set out in the statement of the 

theorem are satisfied. □ 



To describe local substitution rules we prove a theorem that is equivalent to Theorem 2.3. 

Theorem 5.3. Consider a local isomorphism class of substitution tilings given by M and I as 
described in Theorem 5.2, with substitution rules given by s. 

• If I the local substitution rules have s G 

• If I — Ml G the local substitution rules have s E 7? 

• If I + Ml G the local substitution rules have s ^1? + 1 

• For all other values of I, there are no local substitution rules. 

Proof. A substitution rule is local if the label of each tile can be deduced from the geometric 
shape of some bounded patch of tiles around it. 

Section 4.1 shows that the geometric shape of a neighbourhood of the tile can be deduced 
from the projection of its vertices to the window. The window admits a hierarchy of parti- 
tions, with the interval in the partition to which a vertex projects giving the shape of the 
neighbourhood of the tiling near the vertex. Finer partitions give more information about the 
neighbour, in particular giving larger patches of tiling around the vertex. The boundary points 
of the union of all these partitions form the set liyi'^^ + t) U Ily (Z^ + / + t). 

A substitution rule is local if the labelling partition is such that its boundary points precisely 
lie in the set of boundary point for the partition that corresponds to the geometric shape of 
neighbourhoods ny(Z^-|-t)uny(Z^-|-/-|-t). The boundary points are made up of the projection 
of the gifi^s orbits of the points —a and I — b. 

Consider the action of gi^^s on the set Z^ + k, for some point k G M^. The first part of the 
function that takes a point to the strip defined by the subwindow adds a vector in Z^. The 
second part of the function applies a linear map to expand the subwindow to the full window. 
This gives points in the set 

M-\Z'^ + k-s) = Z'^ + M-\k - s). 

Thus £?/,o,s(Z2 + k)cZ^ + M-'^{k - s) 

Now -a G I?, so gifi^s{-a) G Z^ - M'^s. This implies that -M'^s G Z^ U {I? + /). This 
gives two cases s G Z^ and s ^1? — Ml, which we will consider separately. 

If s G Z^ then 

5,,0,.(Z') C Z2, 

thus all elements of the orbit of —a correspond to local labels. Now consider I — b. We have 

gi,o,sil-b) G iP + M'^il-s). As s G iP, there are two possibilities / G and M~^l G iP + l, 
the second is equivalent to I — Ml G Z^. In both cases, the set Z^U(Z^-|-Z) is closed under gi^^s^ 
so the (7i^o,s-orbits of —a and / — b will only give points in U (Z^ + /) and the substitution rule 
is therefore local. If neither of these conditions are satisfied, then the partition will contain 
points which are not in 1? U (Z^ + 1) and the substitution rule will not be local. 

Now consider s G Z^ — Ml, where gi,o,s{~0') G Z^ + /. We must therefore consider the points 
in gi,o,s{'^'^ + 0' to obtain further points on the orbit of —a and the orbit of / — 5. We have 

gi,o,s{^^ + C Z2 + M-^{1 ~s) = I? + M-^l + I. 
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Again this leads us to consider the two cases / G 1? and M~^l + 1^7?. By an analogous 
argument to the one given above, the substitution rule is local if and only if one of these 
conditions holds. 

We reorder the conditions obtained above to emphasize the length of the window (and thus 
the local isomorphism class of the tiling) first, yielding the following cases: 

1. / G and s G U - MZ 

2. I- Ml£l? and s G Z^ 

3. I + Ml el? and seI?-Ml 

These may be further simplified to yield the conditions given in the statement of the 
theorem. In the first case the two sets that can contain s are the same since Z G Z^. In the 
third case we may add I + Ml G I? to the set for s, to obtain s el? + 1 □ 

Finally we show how this result is equivalent to Theorem 2.3. 

Proof of Theorem 2.3. By Theorem 5.3 we have three necessary and sufficient conditions for 
the existence of local substitution rules: 

1. / G Z^ and s G Z^ 

2. I -Ml el? and s G Z^ 

3. 1 + Ml el^ &nd s el? + 1 

The first case holds precisely if both the second and third case hold, so we need only consider 
those. 

If detM = 1, then d- = |ny(s)| and d+ = \Uy{1 - Ml - s)\. If detM = -1 then 
d- = \IIy{s + Ml)\ and d+ = |ny(Z — s)\. There thus remain four cases to be considered, the 
second and third cases from above with det M cither positive or negative. Firstly if det M = 1, 
the second case we have (1 - det{M)/X)\n\,d+,d^ G Uy{1^), as ny(/) = and Uy{MI) = 
det{M)/X\n\. In the third case we have (1 + det{M) / \)\Q\ G ny(Z2), d_ e IIy{1'^ + I) and 

d+ G ny(z -Mi-i? -i) = ny(z2 + z), 

since l + Mlel?. 

Finally consider the case that detM = —1. In this case ny(MZ) has the opposite ori- 
entation of ny(Z). Hence, in the second case it follows that + \MVl\ G ny(Z^), d- G 
Uy{1^+MI) = Uy{1?+1) andd- G Uy{1?+1). In the third case we have |Q|-|Mf^| G ny(Z2), 
d_ G ny(Z2 + l + Ml)= Uy{1^) and d+ G ny(Z2 + 1-1) = Uy{1^). 

We may therefore summarize the cases as in Theorem 2.3: 

• (1 - det(M) /X)\n\,d+,d- e By (Z2) = Z[A] , 

• (l + det(M)/A)|f^|,d+ - - \Q\ G Z[A]. 

□ 
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6 Tilings fixed by a substitution rule 



In this section we address the following question: given a substitution rule, which tilings are 
fixed, ie mapped to (a translated copy of) themselves, by this substitution? Recall that in 
general the predecessor of a substitution tiling is only locally isomorphic to this tiling, and 
often not a (translated) copy. 

Because of corresponding interests in the literature [CMPS93, KvdP96, Yas99], we consider 
not only tilings of the entire line, but also tilings of the positive half-line U {0} (with 
boundary vertex at 0) . 

The results are summarised in the following theorem, which contains the result announced 
in Theorem 2.7. 

Theorem 6.1. Consider a substitution rule a with interval projection structure. Let s,l £ Q'^ 
denote the associated subwindow translation and window vectors. Let X (= det(M) / X) denote 
the contracting eigenvalue of the primitive matrix M G GZ(2, Z), associated to the substitution 
rule. Then 

the number of tilings of the line and half-line that are fixed by the substitution rule are 
as listed in Table 2 and Table 3. These numbers are functions of s, I, M , and the invariant 
window interval (IWI) 

[5l(fl,5lM + n,(i)lcw 

^ A - 1 A - 1 ^ 





dct(M) = 1 


dct(M) = -1 


s U - (M - /)/ 


det(M - /) 


det(M - /) 


s G Z^ U Z^ - (M - 1)1 


det(M - /) + 1 


det(M -I) -I 



Table 2: The number of tilings of the line that are fixed by a substitution rule with an interval 
projection structure. M G GZ(2,Z) is the primitive matrix and s,^ € are the subwindow 
translation vector and window vector associated with the substitution rule. / denotes the 2x2 
identity matrix. 





det(M) = 1 


det(M) = -1 


IWI non-singular 


1 


1 


IWI singular 

on boundary 


1 




IWI singular 
in interior 


2 






Table 3: The number of tilings of the half-line M+ U {0} that are fixed by a substitution rule 
with an interval projection structure with associated primitive matrix M e G/(2,Z), subwin- 
dow translation vector and window vector s,l G Q[A], and invariant window interval (IWI) 
+ Cl^. The IWI is singular if and only if s e (M - 1)1? U (M - I){Z'^ - I). 

Key to the proof of Theorem 6.1 is the observation that the substitution rule induces a 
simple map on the set of windows taking the the window of the original tiling T to the window 
of the tiling aT obtained by application of the substitution rule. 
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Lemma 6.2. Consider a tiling Ti^t with an interval projection structure and window base vector 
t G M^, that admits a substitution rule a as in Theorem 6.1. Then o{'Ti,t) = Ti,Mt-s- 

Proof. We recall that the subwindow of fi;,^ associated with the substitution rule is fi;,r,s := 
\^l^Q + r + s. Then M-'^Tli^r,s = ^i. M-i(r+s) — ^i,t is the window associated to the predecessor. 
As r = Mt — s the substitution rule thus induces the transformation t — > Mt — s. □ 

Finally we prove Theorem 6.1. 

Proof of Theorem 6.1. We first consider tilings of the entire line M. We note that two projection 
tilings are the same up to translation if the projections of the associated window base vectors 
to Y differ by a vector in UyC^^). 

Using the result of Lemma 6.2 we thus have 

UY{t + z) = UY{Mt-s) (2) 

for some z ^1? . 

As before, without loss we consider s G Q'^. Recall that if t € is a window base vector 
for a projection tiling, then any vector in {t + u | u G F} is also a window base vector for the 
same tiling. Since M G Gl{2,Z), we may in fact without loss of generality restrict the solution 
set of (2) to t G and consequently replace this equation by 

t = {M - I)-\s + z), z^l? (3) 

where / denotes the 2x2 identity matrix. 

Viewing (3) to define i as a function of (equivalence classes of) tilings that are invariant 
under this substitution rule are represented by the elements of \t{z) mod 1? \ z ^ 1?}. 

The number of different (closed) window intervals that are invariant under a substitution 
rule is thus equal to 

\{t{z) mod l?\ze = I ((M - = I det(M - 7)|. 

Note that det(M - /) = 1 - T + L», where T = trM, D = dot M and - TA + L> = 0. 

In the argument above we find the number of closed window intervals that are invariant 
under the induced action of a substitution rule. In order to make the relation with tilings we 
need to recall that the window is a half-open interval, which provides the window not only 
with a position and length, but also with an orientation. 

We say that a tiling with interval projection structure (and its window) is singular if 
dni^t + F n Z2 / (or equivalents if dQ.i^t n ny(Z2) ^ 0) and non-singular otherwise. 

The base vector t of a singular window satisfies 

ny(i) GnY(z2)uny(z2)-nY(/). (4) 

Since I\-y{1) G Q[A] we can without loss choose / G and subsequently solve for i G in the 
case of singular tilings, rewriting (4) as 

t el? yji? - I. 

It thus follows that within a local isomorphism class there are one {I G 1?) or two (I 1?) 
different tilings (up to translation) that are singular. 
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Consequently, there exists a singular tiling that is candidate to be invariant {at e l?\J7? — l 
that satisfies (3)) if and only if 

s - {M - 1)1. 

If s ^ 1? VM? — (M — T)l there exists no singular invariant tiling and the orientation of 
the window needs no attention as the boundary of the window of an invariant tiling cannot 
intersect the lattice 1? . For such substitution rules the number of different invariant tilings is 
thus exactly equal to | dct(Af — /)|. 

In the singular case it is important to compare the orientation of the window before and 
after substitution. Note that for each s ^ 1? VM? — [M — 1)1 there is precisely one choice 
of t (mod Z^) satisfying (4) and thus precisely one singular tiling that is a candidate to be 
invariant under the substitution rule. 

If det M = 1 the windows before and after the geometric transformation associated with 
substitution have the same orientation. For each singular window base vector t, we thus find 
two different singular invariant tilings: one for each choice of the side of closure of the window. 
Consequently, in this case the number of invariant tilings is | det(M — /)| + 1. When det M = —1 
the orientation of the singular tiling changes so the substitution rule cannot fix it and the total 
number of invariant tilings is | det(M — /)| — 1. 

The results are summarized in Table 2. 

It now remains to analyse the number of fixed tilings of the half- line R"*" U {0}. In this case 
there is no way to allow for translations, so that invar iance implies 

t = Mt-s ^ t= {M - I)~^s, 

with s G Q^. The corresponding invariant window interval (IWI) is [^x^! + ^1 C W . 
We first note that Q lies inside the IWI. In fact, is the unique vertex that is also a 

vertex for all the expansion predecessors. Hence must be the boundary vertex. 

If the window translation vector s of a substitution rule satisfies s (M — 1)1? U (M — 
— I) then the IWI is non-singular and there is a unique invariant tiling of the half-line 
M+ U {0}. 

If the window is singular then we must consider the sign of determinant of the matrix M 
and the location of (projections) of the singular lattice point(s). If the tiling admits a vertex 
that is the projection of a singular lattice point (lying on the boundary oiV + Q^t) then we 
call this vertex singular. 

If the IWI is singular, but the tiling does not contain a singular vertex (as all singular 
lattice points project to the negative part of the half- line), there is no interference with the 
tiling and thus there is one unique invariant tiling of the half-line U {0}. 

If is at the boundary of the IWI it is a singular vertex. As serves as the boundary point 
we have no choice with regard to the side on which to close the window. Moreover it turns 
out that can be on the boundary only if det(M) = Consequently, there is one unique 
invariant tiling of the half-line U {0} . 

If the tiling for the IWI has a singular vertex that does not lie on the boundary we have 
to distinguish two cases, depending on the sign of det(M). If the sign is positive there are two 
invariant tilings, corresponding to the two choices of the side where the window is closed. If 
the sign is negative there are no invariant tilings. 

The results are summarized in Table 3. □ 
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7 Interval exchange transformations 



Our results on the existence of substitution for one-dimensional tilings with an interval pro- 
jection structure are intimately related to renormalization properties of two and thiee- Interval 
Exchange Transformations (lETs) . In this section we discuss this relationship and reformulate 
our main result in terms of lETs. 

lETs are defined as follows. Let / C M be a bounded half open interval that is partitioned 
into a set of half open subintervals {/q | a G .4} indexed by some alphabet .4 of d > 2 symbols. 
An interval exchange transformations is a bijection from I to I which is a constant translation 
on each of the subintervals la, thus inducing a permutation of the subintervals covering /. 
An interval exchange transformation / is determined by combinatorial and metrical data as 
follows: 

(1) A pair tt = (ttojTTi) of bijections tt^ : v4 — > {1, . . . ,d}, with e G {0, 1}, describing the 
ordering of the subintervals la before and after the transformation is applied. One may 
represent this permutation as 



(2) A vector /x = {iJ,a)aeA with positive entries, where /Xq, denotes the length of the subinterval 



We call p = TTi o ttq the monodromy invariant of vr. It should be noted that we can normalize 
the combinatorial data by choosing ^ = {1, 2, . . . , d} and ttq = id so that tti = p and tt = {id,p). 
We denote an interval exchange transformation (lET) / by the pair / = (tt, yu). 

The study of interval exchange transformations has a rich history, mainly motivated by 
the study of translation surfaces. For more details, we recommend the recent survey by Viana 
[Via05]. 

lETs appear in the study of one-dimensional tilings with an interval projection structure, 
via the map Ily o whose constituents have been introduced in Section 3.2. Orbits of this 
map represent projections of subsequent vertices of the staircase to the window Q.i^t C Y . From 
Definition 3.4 it follows that Ily o ^ is an lET on two or three intervals covering the window 
Such lETs are considered in detail in [FHZOl, FHZ03, FHZ]. 

The type of lETS arising in this context of Section 3.2 are precisely the irreducible two 
and three-IETs. An lET / on an interval / is called irreducible if there does not exist any 
subinterval I of I for which f\j is an lET on /. These lETs may be represented by the following 
two-parameter family of maps. Let / = [0, 1) denote the half open unit interval and a,b E I 
two points such that < a < 6 < 1. Then we define 



The map fa,b is a two-IET if a = 6 and a three-IET otherwise. 

In the study of lETs, considerable attention has been devoted to the return map Rjf that 
an lET / = (tt, /j,) induces on a half open subinterval I C I. It turns out that the return map 




where aj = tt^ ^{j) and j G {1, 2, . . . , 



d}. 



la. 



' X + 1 — a if X G [0, a), 
fa,bix) = < x + l-a-b ifxG[a,6), 
X — b if X G [6, 1), 
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is again an lET, with permutation 7f of a partition Jl oi I that is imphcitly defined by the 
dynamics of /. In order to compare the lET f on I with the return map that it induces on 
/ C /, it is natural to identify / with / (by means of a uniform expansion with factor |/|/|/|, 
an appropriate translation, and - in case the orientations of / and / are opposite - a reflection). 
We refer to this mapping between lETs as the Rauzy-Veech renormalization induced by /, 
after [Rau79, Vee78]: 

Rj{7r,fj,) = {W,-fl). 

More information about the dynamics of the original lET is retained by considering an 
extension of the Rauzy-Veech renormalization that includes suspension data. In addition to 
considering the induced lET on the subinterval /, one may choose to keep track of the paths 
that suspended orbits take in 7, while being taken by iteration of / from I through I\I before 
returning to /. The suspension data may be represented by an element of the monoid A* 
consisting of all finite words in the alphabet A that is used to label the subintervals of /. 
We thus define the Rauzy-Veech renormalization on the lET including suspension data as 

'Rj{iT,n,v) = {W,JI,v), 

where the suspension data are represented hy v = {vi, . . . ,Vn) with Vi an element of the monoid 
A*. 

Recall that a morphism, or substitution rule, a : A* ^ A* associates to each element of 
the alphabet A a finite word over A and acts on words in w{ai,a2, ■ ■ ■ , (Xd) E A* as 

aw{ai,a2, ...,ad)= w{a{ai),a{a2), cr{ad)). 

We address the question for which lET (tt, /x) and subinterval / there exists a refinement jl of 
the partition on which vr induces a permutation tt, so that the Rauzy-Veech renormalization 
mapping fixes (rr,/!) and induces a morphism a on the corresponding suspension data v = 
{vi,... , vg) with Vi G A* where A is an alphabet with d element: 

nj{TT,fi,v) = {TT,fl,av). (5) 

Theorem 7.1. Let fa^b = /^) a two or three-interval exchange transformation on the 
unit interval I with discontinuity points < a < 6 < 1. Then there exists a refinement jl of fi, 
with permutation tt trivially induced from n, such that the Rauzy- Veech renormalization with 
respect to a subinterval I induces a morphism a on the suspension data v, as in (5), if and 
only if there exists a quadratic unit X, such that 

• |A| < 1, 

• its algebraic conjugate is positive, 

• A(Z-FZ(a- l)/6) =Z-hZ(o- l)/6, 
andl = t + \I d I with t G Q[A]. 

Proof. This theorem is a consequence of Theorem 2.1 and Theorem 2.2. Prom the discussion in 
Section 4 and Section 5 it follows that the existence of a substitution rule for a one-dimensional 
tiling with interval projection structure is equivalent to the existence of an induced morphism 
on the suspension date of Rauzy-Veech renormalization for a two or three-IET. The subinterval 
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/ corresponds to the subwindow associated to the predecessor tihng and the refinement /x to 
the labelhng partition of the window, and the morphism a to the substitution rule. 

Let us consider an lET with interval length 7 > and directed one-dimensional vectors 
a < < /3, satisfying \a\ + \ f3\ > 7, as considered in Section 4. This corresponds to an lET on 
the unit interval / with discontinuity points < o < 6 < 1, where a = (7 + a)/7 and h = jd/^y. 

In the setting of Section 4, a and (3 are considered as projections (to Y) of generators of 
a two-dimensional lattice isomorphic to I?. The projection of this lattice to Y is aL + (5Z. If 
there exists a lattice automorphism that acts as a uniform contraction by a factor A on y and 
as a uniform expansion by a positive factor containing the interval, it follows that A must be a 
quadratic unit, and the point set a'Z + (37. must be invariant under multiplication by A, which 
implies that 

ny(z2)(Z6 + TL{a - 1)) = Z6 Z(a - 1). 

This implies in particular that there exist integers x, y such that \ = x ^ y{a — \) jh and hence 
A G Z -I- Z(a - 1) /6, but also that A = x6/(a - 1) -I- y so that (a - 1)/6A G Z -|- Z(a - 1) jh which 
in turn implies (as A Z) that (a — must be a quadratic algebraic number. 

On the window and subwindow / we have the following conditions from the corresponding 
conditions for substitution rules in Theorem 2.1 and Theorem 2.2. First, the length of the 
window interval 7 must be in Q[A], and second, so must the distances between any of the 
boundaries of the subwindow interval to any boundary of the orginal window interval. In 
terms of the normalized lET this implies that the distances between boundary points of / to 
any of the boundaries of / must be in Q[A], i.e. 

dl c Q[A] n /. 

□ 

Note that the conditions on A are necessary and sufficient conditions for the existence of 
a subinterval / on which the Rauzy-Veech renormalization induces a substitution rule on the 
suspension data, in correspondence with Theorem 2.1. The final condition corresponds to the 
one set out in Theorem 2.2 characterizing the countably infinite number of appropriate choices 
for the subinterval /. 

Theorem 7.1 relates to some results on renormalization by [BC97], where Rauzy-Veech 
renormalization on specific subintervals is considered, without considering the corresponding 
action on suspension data. 

We finally illustrate the Rauzy-Veech renormalization on the Fibonacci lET. 

Example 7.2 (Rauzy-Veech renormalization of the Fibonacci lET). To illustrate the 
connection between renormalization of lETs and substitution rules for one-dimensional tilings 
with an interval projection structure, we consider the Fibonacci lET, whose orbits correspond 
to Fibonacci tilings. We focus in particular on the derivation of the non-local substitution rule 
for the Fibonacci tilings, presented in Example 2.5. 

The Fibonacci two-IET on / = [0,1) is defined as jaf) with a = b = . Choosing 
the quadratic unit A = —a, the conditions on A of Theorem 7.1 are satisfied. We choose 
I = XI + t as the central subinterval of the unit interval /, with A = — a quadratic unit 

and t = ^'^/^ G Q[A] \ Z[A]. Satisfying hereby the final condition of Theorem 7.1, we thus 
expect to have some refinement // of the partition /i for which the Rauzy-Veech renormalization 
to / induces a morphism (substitution rule) to the suspension data. 
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ai 





# 



a2 



ai 



Figure 9: Illustration of the iterative construction of a refined partition of the interval / for which 
the Rauzy-Veech renormalization on the central subinterval / of length = — A induces 

a substitution rule on the suspension data for the Fibonacci lET, defined in Example 7.2. 
The suspension data is represented as a graph over the partition of the interval. Because 
the contraction factor A < 0, in comparing the lET on / to the induced lET on /, we need 
not only rescale and translate, but also perform a reflection in the line. We start with the 
Fibonacci lET. schematically indicated in the left graph under (i), and the induced lET on / 
indicated in the right graph. It follows that the Rauzy-Veech renormalization fixes the lET 
i?j(7r, /u) = (tt, /i), but it does not induce a morphism on the suspension data, since in the 
graph over the longest subinterval of / one distinguishes two different graphs: ab and ba, 
corresonding to different suspensions. In order for the renormalization to induce a morphism 
on the suspension data, corresponding to each marked subinterval of / there should be unique 
suspension data. We thus need at least make a refinement of the a subinterval of / into two 
subintervals which we label ai and 02 (left of (ii)). If we examine the tower on the induced 
partition {la^ , I-^^Ij;} of /, we observe that the graphs over and are well defined, but 
that there is a problem with the suspension over Jj-. We subsequently consider the additonal 
refinement depicted at the left of (iii) , where the subinterval is partitioned in two subintervals 
/foj and Ib2- With this partition, we observe at the right figure under (iii), that the suspension 
data on the induced lET on / now is unique on each labeled subinterval of /. From the graph 
we read the morphism a{ai, a2, 61, 62) — (&iai, a2^^2, a2, ai). This is precisely the substition rule 
mentioned in Example 7.2. 
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In Figure 9 we illustrate how the refinement of the partition can be constructed by ex- 
amining the return map to I, including suspension data. We note that due to the fact that 
A < 0, the orientation of the half open interval / is opposite to that of /. This implies that 
in the identification of the lET induced on I by the return map with I, we need to take into 
consideration a refiection. As the boundary points of / lie in Q[A] \ Z[A], the substitution rule 
is non-local, in the sense of the discussion in Section 2. 
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